Microwave conductivity of d-wave superconductors with extended impurities 



Tamara S. Nunner and P.J. Hirschfeld 
Department of Physics, University of Florida, Gainesville, FL 32611-8440, USA 

(Dated: February 2, 2008) 

We investigate the influence of extended scatterers on the finite temperature and finite frequency 
microwave conductivity of d-wave superconductors. For this purpose we generalize a previous treat- 
ment by Durst and Lee [Ref. 22], which is based on a nodal approximation of the quasiparticle 
excitations and scattering processes, and apply it to the analysis of experimental spectra of YBCO- 
123 and BSCCO-2212. For YBCO, we find that accounting for a slight spatial extension of the 
strong scattering in-plane defects improves the fit of the low temperature microwave conductivity 
to experiment. With respect to BSCCO we conclude that it is necessary to include a large con- 
centration of weak-to-intermediate strength extended scatterers, which we attribute to the out-of 
plane disorder introduced by doping. These findings for BSCCO are consistent with similar analyses 
of the normal state ARPES spectra and of STM spectra in the superconducting state, where an 
enhanced forward scattering has been inferred as well. 
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I. INTRODUCTION 

Early in the debate over the symmetry of the su- 
perconducting order in the cuprates, a rather convinc- 
ing picture of the microwave properties of the YBCO- 
123 system was put forward by Bonn et al.^ and later 
placed on a microscopic foundations 3 ^. Crucial to this 
interpretation is the observed collapse of the d-wave 
nodal quasiparticle scattering rate as the system becomes 
superconducting£*£*£, leading to a dramatic rise in the 
conductivity with decreasing temperature. This rise is 
cut off when the inelastic mean free path becomes com- 
parable to the elastic one, and the conductivity subse- 
quently decreases because of the vanishing nodal carrier 
density as temperature tends to zero. One consequence of 
this picture is that the resulting conductivity peak should 
be suppressed and occur at higher temperatures in dirtier 
samples. In addition, the conductivity should approach 
the universal value cro = e 2 vp/hv2 for zero temperature 
and zero frequency as predicted by P. A. Lee& for the case 
of isotropic scatterers, where Vp is the Fermi velocity and 
v-i the gap velocity at the node. 

If one extracts Vp / u 2 from thermal conductivity^ or 
angle resolved photoemission (ARPES) measurements^, 
one finds that the universal value for both BSCCO and 
YBCO crystals should be about a = 5x 10 5 O _1 m _1 . In 
YBCO, the residual value of the conductivity for lu,T — >0 
is difficult to determine, but appears to be approaching 
3-4 x Co in the best crystals^. The peak in the conduc- 
tivity occurs around 25K with an amplitude of approx- 
imately 100 x (To for the lowest frequency measured. In 
BSCCO, the peak is located around 20K, but is only 
about 20% higher than the apparent residual valued of 
8-IOxctq- Virtually no frequency dependence is seen in 
the measured microwave frequency ranged, suggesting a 
very dirty material, in apparent contradiction-within the 
"standard" scenario-with the low-temperature peak po- 
sition. The longstanding puzzle of the low temperature 
microwave peak together with indications of dirty limit 



behavior have been analyzed as evidence for absorption 
into a collective mode off resonance at low frequencies 13 , 
as well as a consequence of nanoscale inhomogeneityi^. 

In this work we argue that the temperature depen- 
dence of the conductivity can be more naturally under- 
stood in terms of the effect of extended scatterers present 
in the BSCCO crystal. Current generation crystals are 
made typically with excess Bi, deficiencies of Sr and Ca, 
and variable O content; cation substitution is thought 
to occur frequently. Some aspects of this defect distri- 
bution have been discussed recently by Eisaki et ali^. 
The net result of these defects is not only to dope the 
nominally stoichiometric BSCCO crystal (pure BSCCO 
would be an insulator), but to provide a relatively slowly 
varying potential landscape for quasiparticles moving in 
the Cu02 planes. The effect of these extended scatter- 
ers with respect to the normal state has recently been 
discussed by Abrahams and Varmai& and it has been 
pointed out by Zhu et tdm± that the broad momentum 
space peaks observed in Fourier transform STM stud- 
ies of BSCCOiS*i2iSli can only be explained by poten- 
tial scatterers with finite range. In a further application 
of this notion to ARPES, Zhu et al 21 showed that a 
large concentration of impurities with potentials peaked 
in the forward direction could be present without sub- 
stantially broadening quasiparticle states except near the 
node. Since the microwave conductivity is dominated by 
nodal quasiparticles, it is clearly important to ask what 
the effects of extended or forward scatterers are in this 
case. 

Since the work of Durst and Lee^ we know that the 
residual conductivity in the presence of extended scat- 
terers can be much larger than the "universal" value <To. 
This might account for the large value of the microwave 
conductivity observed in the BSCCO-2212 system at low 
temperatures. To make this case, however, one needs 
to examine the influence of a finite scattering range at 
finite temperatures and frequencies. We have therefore 
generalized the analysis of Durst and Lee in this way and 
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applied this treatment to the analysis of experiments on 
YBCO and BSCCO. With respect to YBCO we find that 
consideration of slightly extended strong scatterers pro- 
vides a better fit to the low temperature microwave con- 
ductivity than pointlike strong scatterers. For BSCCO 
we conclude that it is necessary to include a large con- 
centration of weak-to- intermediate strength extended im- 
purities in addition to the strong in-plane defects which 
are responsible for the unitary scattering resonances ob- 
served by STMM. 

The outline of the paper is at follows. In section [B] we 
describe the model and derive expressions for the self- 
energy and vertex function for extended scatterers. Our 
approach, which is based on an extension of the work by 
Durst and LeeSS, aims at treating scattering potentials 
with an extension of a few lattice spacings at maximum 
and is therefore in the opposite limit from semiclassical 
calculations where the impurity potentials extend over 
a few coherence lengthsS^S^. In section IIIII we apply 
our treatment to the microwave conductivity of YBCO. 
We show that the consideration of slightly extended in- 
stead of pointlike strong potential scatterers improves the 
agreement with the experimental spectra. In section ITvl 
we address the microwave conductivity of BSSCO and 
demonstrate that it is necessary to include a large con- 
centration of weak extended scatterers in order to explain 
the experimental spectra. A good fit is obtained based 
on a realistic disorder model for BSCCO which contains 
weak extended scatterers in addition to strong pointlikc 
impurities. Finally, in section IV, we present our conclu- 
sions. 



II. TREATMENT OF EXTENDED 
SCATTERERS 

For low temperatures and low frequencies the quasi- 
particle dispersion of a d-wave superconductor can be 
linearized around the nodes. The resulting quasipar- 
ticle spectrum has the form of a Dirac cone, whose 
anisotropy is determined by the ratio Vfjv% of the Fermi 
velocity Vf = de k /dk = 2\f2t and the gap velocity 
v 2 = dAk/dk = Ao/v2, where t is the nearest neighbor 
hopping parameter and Ao is the maximum gap value 
and we have set a = % = 1. At low temperatures and fre- 
quencies, quasiparticle excitations are restricted to small 
regions around the nodes. Therefore momentum trans- 
fer between quasiparticles is limited to four wavevectors 
which connect the four nodes and include intranode and 
internode scattering processes (see Fig.^l. Consequently 
a momentum dependent impurity potential V{k) can be 
represented by three parameters V\ , V 2 and V3 which cor- 
respond to the respective momentum transfers. Based 
on these approximations and treating impurity scatter- 
ing in T-matrix approximation, an expression for the 
microwave conductivity has been derived by Durst and 
Lee22. They found that vertex corrections, which arise 
from the momentum dependence of the impurity poten- 



tial, induce a dependence of the zero-temperature and 
zero- frequency limit of the conductivity on the impurity 
potential and the impurity concentration. Contrary to 
the case of pointlike scatterers, no universal value of the 
electrical conductivity exists therefore in case of extended 
scatterers. Durst and Lee, however, did not further ex- 
plore the frequency and temperature dependence of the 
microwave conductivity. Here, we slightly modify their 
approach and evaluate the conductivity for finite frequen- 
cies and temperatures. Furthermore we consider the ef- 
fect of combining different types of scatterers and calcu- 
late the microwave conductivity for BSCCO based on a 
realistic disorder model. 



A. Self-energy 

Before proceeding to two-particle quantities like the 
microwave conductivity, it is instructive to focus first 
on single-particle properties like the single-particle self- 
energy. Using the Nambu notation, the single-particle 
self energy in a superconductor can be decomposed as: 



£(fc,w) = ^S Q (fc,Lj)f Q , 



(1) 



where f a are the Pauli matrices and fo is the unit matrix. 
Treating impurity-scattering in T-matrix approximation 
gives rise to the following self-energy 



T,(k,uj) = riiT kk (uj), 



(2) 



where rii is the impurity concentration and Tkk(to) is the 
diagonal element of the T-matrix 

Tifcfc'M = Vkk'h + Vkk"f$G(k", u))fw> k ' M ■ (3) 



The self-energy u>) has to be solved self-consistently 
in combination with the single-particle Green's function 

Gik,^)- 1 = G {k,uj)- 1 -t(k,uj) (4) 

where the unperturbed Green's function is given as 

cjtq + A fe fi + e fe f 3 



Go(fc,w) 



(5) 



Following the approach of Durst and Lee2£, we reduce 
the impurity scattering potential to the four wave- vectors 
connecting the nodes, i.e., Vkk' is replaced by a 4 x 4- 
matrix in nodal space 
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where V\ , V 2 and V3 are the values of the impurity poten- 
tial at the wavevectors connecting the nodes, see Fig. ^ 
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FIG. 1: For low temperatures and low frequencies the mo- 
mentum transfer between quasiparticles is basically limited 
to the four wavevectors connecting the nodes of a d-wave su- 
perconductor. Therefore the momentum dependent impurity 
potential V(k) can be represented by the values at the respec- 
tive wavevectors, i.e., by three parameters Vi, V2 and V3. 



Using this simplification, the impurity potential can 
be pulled out of the integral and the T-matrix becomes 
a 4 x 4-matrix in nodal space 

Tjj> (w) = Vjfh + lG(u)h v n"Tj"j> M, (7) 

3" 

where Ig(^>) is the integral of the single-particle Green's 
function over one quarter of the Brillouin zone 



d 2 k , 



I G {uj)f . (8) 



Jo (2^) 2 

This allows for an analytical solution?- of the T-matrix 
^•'=^3+7^0 (9) 

with 

V \ mn ( -I G (UJ)V 2 
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T "' \1-Ig(u) 2 V 2 



ii' v x ^r; i / jj> 

(10) 

where the denominators have to be calculated as inverse 
matrices in nodal space. This gives for the Eo-component 
of the self-energy: 



E (w) 



4/gM 
1 



l-I&CwXVx-Va) 2 



(11) 



l-J G (a/)(Vi - 2V 2 + V 3 ) 2 1- J*(w)(Ki + 2V 2 + v 3 y 

For an isotropic impurity potential, i.e., Vi=V2=V3=V 
this expression simplifies to: 



4/ G (w) 



1 



1 



1 - 16/ G H^ 



(12) 



which implies that the self-energy for a momentum- 
dependent impurity potential in Eq. <| 1 1 f) can be decom- 
posed into a sum of self-energies corresponding to three 
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FIG. 2: Imaginary part of the self-consistently calculated self- 
energy Eo(w) for an anisotropic impurity potential character- 
ized by the three parameters Vi/t = 100, V^/i = 66 and 
Vz/t — 52 (solid line). The positions of the resonances co- 
incide with the resonances of the self-energies for isotropic 
impurity potentials V/t — 71, V/t = 12 and V/t = 5. Here 
Ao/t = .29 and m = .00002 have been used. 



different isotropic impurity potentials 
E (w) = TgS°(j(Vi + 2V 2 + V 3 ),u) 



(13) 



+4 s °(J(Vi - 2V 2 + V 3 ),u) + 2£ni(V 1 - V a ),u>) . 

Consequently, the self-energy for an anisotropic impu- 
rity potential in the nodal approximation contains up to 
three resonances corresponding to the different impurity 
strengths V=(Vi + 2V 2 + V 3 )/4, V = (Vi - 2V 2 + V 3 )/4 
and V ={V X - V 3 )/4, see Fig. El 

If the scattering strength of the impurities is weak, 
they can be treated within Born approximation and the 
self-energy of an extended weak impurity becomes: 



£ (w) 



-n l {V 2 + 2V 2 + V 2 )I G {u), (14) 



i.e., the self-energy for an anisotropic impurity potential 
is identical to the sel f-energy for an isotropic impurity 
potential with V = yJV? + 2V 2 2 



V 2 
v 3 ■ 



B. Microwave conductivity 

In linear response the electrical conductivity is given 

Imn ret (0,T) 



<r(n,T) 



(15) 



where II rc t(f2,T) is the retarded current-current correla- 
tion function, which can be obtained by analytical con- 
tinuation from 



2 2 
ev f 



U(ifl) = ^ Tr[G(k, iu)G(k, iu+iQ)f(k, iuj + ifl)} , 

1 ■ 1 

lUJ,k 

(16) 
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where T(fc, iu + in) is the vertex function, which for a 
d-wave superconductor arises entirely from the momen- 
tum dependence of the impurity potential and will be 
calculated here as the sum of ladder diagrams. 

For simplicity we neglect the £3-component of the self- 
energy and keep only the So-component. This approxi- 
mation becomes exact in the unitary limit for isotropic 
scattering V\ = V% = V3 and for purely forward scatter- 
ing ^2 = ^3 = because the £ 3 -component vanishes in 
theses limits and also in the Born approximation for all 
scattering potentials. In terms of diagrams this means 
that all summations are reduced to diagrams with even 
number of impurity lines. In order to obtain a conserving 
approximation^ for the conductivity, the vertex correc- 
tions also have to be restricted to diagrams containing 
even numbers of impurity lines, see Fig. |31 This is a 
modification to the treatment by Durst and Lee22, who 
have summed all ladder diagrams both with even and odd 
number of impurity lines. This modification is necessary 
when E 3 is neglected because otherwise the vertex correc- 
tions violate analyticity at finite to and T. Furthermore, 
this modification leads to a simplification to the expres- 
sion of the conductivity obtained by Durst and Lee2£, 
see below. Note, however, that our treatment still agrees 
with Durst and Lee2£ in the Born approximation, which 
is not affected by this modification and also in the T- 
matrix approximation for the limit ui —* and T — > 0, 
i.e., the limit Durst and Lee have focused on. Differences 
arise only for finite temperatures and frequencies in the 
T-matrix approximation. 

• A *"* 



where lB(ioj,i£i) is the momentum integrated particle- 
hole bubble 



Igiiw, in) 



d 2 k - 



W 



G{k,iu))G{h,iw + iti). (19) 



Note, that here we calculate all integrals Ig(iu},iSl) and 
Ig (w) numerically instead of replacing the elliptical inte- 
gration area by a circle as in Ref. 2^, because this approx- 
imation becomes inaccurate for small values of the gap, 
i.e., close to T c , where we parameterize the temperature 
dependence of the gap in the following wayi22 



A (T) = A tanh(aV?VT-l) 

with a = 3.0. 

The vertex function is given by 



(20) 



(21) 



{T^{iLu)T^(iuj + in) + T^(iu)T^{iLU + m) 
~T? 3 (ioj)T? 3 (iuj + in) - T? 3 (iw)T? 3 (ioj + in)) . 

After analytical continuation the microwave conductivity 
can be expressed as^ 



a(n) 
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e v f 



riF^Lo) — uf{oj + n) 



(22) 



2tt 2 v 2 J n 

(Re J(u> - iS,u + n + i8) - Re J(lj + id, w + fi + iS)) . 

In Born approximation one arrives at a very similar 
expression for the conductivity, but the vertex function 
7(iw, in) in Eq. (|18|) is replaced by the simpler from 




FIG. 3: T-matrix approximation for the self-energy £0 and 
the microwave conductivity a. Because we neglect the fa- 
component of the self-energy, which corresponds to neglecting 
all diagrams with odd number of impurity lines, all diagrams 
with odd number of impurity lines have to be excluded from 
the vertex corrections as well. 



Summing up all ladder diagrams with even number of 
impurity lines one arrives at the following expression for 
the current-current correlation function: 



2 

n(zf2) = t^L y jfiw.ifi) 
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with 
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Ib{iijJ, in) 



1 — 7(20;, in)Is{iuj, in) 



(17) 



(18) 
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So far we have focused on the effect of impurity scatter- 
ing, which is appropriate for low temperatures and low 
frequencies. At higher temperatures, however, it is es- 
sential to take into account inelastic scattering processes 
as well, like e.g. quasiparticle quasiparticle scattering or 
scattering off spin fluctuations. These inelastic scatter- 
ing processes are suppressed below T c due to the opening 
of the superconducting gap in the quasiparticle excita- 
tion spectrum and therefore the contribution of inelas- 
tic scattering increases rapidly when T c is approached 
from the low temperature side. A full treatment of in- 
elastic scattering is beyond the scope of this paper. It 
has been pointed out by Walker and Smith 2 — that the 
contribution of quasiparticle quasiparticle scattering to 
the transport lifetime is exponentially suppressed for low 
temperatures because only Umklapp scattering processes 
can decay the current and a finite excitation energy Ajj 
is necessary to permit an Umklapp scattering process for 
a realistic Fermi surface. Thus, we include the effect of 
inelastic scattering by simply adding the inverse trans- 
port lifetime t~*j(T), which has been obtained by Duffy 
et alm^ by extracting the Umklapp scattering processes 
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from scattering of quasiparticles off spin fluctuations, to 
the imaginary part of the self-energy T, (oj) in Eq. Hll|) 
or Eq. |[Tl|l 

EH = SoM- l (2r incl (T))- 1 (24) 

Note, that our simplified treatment of t-~^ completely 
neglects the frequency dependence of inelastic scattering 
and therefore limits our approach to small frequencies in 
the microwave regime, and prevents us from calculating 
the conductivity in the THz-range. 



III. COMPARISON WITH EXPERIMENTAL 
SPECTRA OF YBCO 

The microwave conductivity has been investigated in 
detail for pointlike scatterers within the self-consistent 
T-matrix approximation 2 - and good agreement with the 
experimental data of YBCO has been found. The tem- 
perature dependence of the microwave conductivity for 
pointlike scatterers, see also Fig. 0] can be summarized 
in the following way. For zero temperature and zero fre- 
quency the conductivity approaches a universal value 8 
(To = e 2 Vf/hv2 due to the fact that at zero tempera- 
ture impurities give rise to a finite quasiparticle density 
of states while at the same time they reduce the quasi- 
particle lifetime. At low temperatures the conductivity 
rises with increasing T due to an increase in the number 
of excited quasiparticles. The exact temperature depen- 
dence is determined by the density of states in a d-wave 
superconductor and the frequency dependence of the im- 
purity scattering rate. Starting from the opposite side, 
i.e., decreasing the temperature below T c the conductiv- 
ity also increases rapidly because inelastic scattering is 
suppressed below T c due to the opening of the super- 
conducting gap in the quasiparticle excitation spectrum. 
This leads to the formation of a peak at intermediate 
temperatures, whose position is determined by the mi- 
crowave frequency and the impurity scattering strength. 
This peak moves to higher temperatures and its ampli- 
tude decreases with increasing microwave frequency and 
impurity concentration. 

Our best fit to the experimental spectra of YBCO 11 
using pointlike strong scatterers is displayed in Fig. ^ 
Commonly used parameters for YBCO are: Ao = 400A" 
for the gap maximum, v//i>2 = 14 for the anisotropy of 
the Dirac cone and T c — 88.7A 911 . In order to com- 
pare our theoretical curves to the experimental data the 
value of the universal conductivity do — e 2 / '(hTT 2 )v / / 'v-i 
has to be translated into a three dimensional conductiv- 
ity which can be done via aQ D = aQ D n c , where n c is the 
number of Cu02-planes per unit length in the c-direction, 
which is n c =1/(5. 9A) for YBCO. Because a ~ A~ 3 the 
conductivity a is very sensitive to the value of the pen- 
etration depth A, which has recently been measured 30 
as considerably smaller than previously published in the 
literature 3 ^, i.e., A = 1030 ± 80A instead of A ~ 1550A. 
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FIG. 4: Fit to the experimental spectra of YBCO (reproduced 
from Ref. I 111) using pointlike strong scatterers with V = WOt 
and a concentration of rn — .000035. The magnitude of the 
conductivity has been scaled by a factor of .42, which corre- 
sponds to assuming a penetration depth of 1200 A. 



This would increase the published values^ of the mi- 
crowave conductivity by a factor of four. Indeed, it turns 
out that we obtain the best fit to the microwave conduc- 
tivity of YBCO when we assume the absolute values of 
the conductivity to be roughly twice the previously pub- 
lished data 11 , which would correspond to a penetration 
depth of approx. 1200A. Therefore we allow ourselves 
the freedom to scale our calculated curves by roughly a 
factor of 1/2 when comparing to the experimental pub- 
lished values (exact scaling factor is stated in the figure 
captions). 

As can be seen from Fig. 0] the assumption of point- 
like scatterers can reproduce the temperature and fre- 
quency dependence of the microwave conductivity of 
YBCO quite well (see Refs. 000). The largest discrep- 
ancy arises for low temperatures and low frequencies, 
where experimentally a nearly linear increase of the con- 
ductivity with temperature is found whereas the theory 
based on pointlike scatterers predicts a quadratic tem- 
perature dependence. It has been suggested by Hettler 
and Hirschfeld 32 that the theoretical lineshape becomes 
more linear at low temperatures when the suppression of 
the order parameter surrounding a strong pointlike scat- 
terer is taken into account. This low temperature be- 
havior has been attributed to the formation of a second 
resonance in the self-energy So at low frequencies. It is 
intriguing to note that we find a similar resonance in the 
self-energy for slightly extended strong potential scatter- 
ers, see Fig. |21 and therefore it is interesting to investigate 
whether the presence of slightly extended potential scat- 
terers can also explain the linear T-dependence of the 
microwave conductivity for low temperatures. Our best 
fit to the experimental spectra of YBCO using slightly 
extended potential scatterers is displayed in Fig. [SJ Ob- 
viously the consideration of extended scatterers consider- 
ably improves the agreement with the experimental data 
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FIG. 5: Fit to the experimental spectra of YBCO (repro- 
duced from Ref. flit) using slightly extended strong scatterers 
with Vi = lOOt, V2 = 92t, Vs = 82t and a concentration of 
rii — .000015. The magnitude of the conductivity has been 
scaled by a factor of .425 which corresponds to assuming a 
penetration depth of 1200A. 



at low temperatures. 

Surprisingly, the concentration of extended scatterers 
used for the fit in Fig. [5] is even lower than the concen- 
tration of pointlike scatterers we have used for the fit in 
Fig. [J] Generally one would assume that due to the for- 
ward scattering character of extended impurities a larger 
concentration is necessary to obtain a similar scattering 
rate as for pointlike impurities. To gain more insight 
into this unusual behavior we show in Fig. the mi- 
crowave conductivity for different extensions of the scat- 
tering potential at two of the experimentally measured 
frequencies, i.e., 1.14 GHz and 13.4 GHz. Increasing 
the forward scattering character of the impurity potential 
slightly from V x = V 2 = V 3 = 100* to V x = 100*, V 2 = 90i 
and V3 = 80£ essentially reduces the height of the peak in 
the conductivity at 1.14 GHz and makes the low temper- 
ature increase more linear. For this small deviation from 
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FIG. 6: Microwave conductivity for strong scatterers with 
varying degree of forward scattering. Impurity concentration: 
m — .000035, maximum gap: Ao = .29i. Left panel: iv = 
1.14GHz, right panel: u> — 13.4GHz. Inset: Self-energy for 
the respective scattering potentials. 



isotropic scattering, vertex corrections are small and the 
variation of the conductivity can be attributed to the 
formation of a second resonance in the self-energy at in- 
termediate frequencies, see inset in Fig. (see also dis- 
cussion in Sec. Ill A|) . The nearly linear increase of the 
self-energy below the second resonance causes the more 
linear T-dependence of the conductivity at low tempera- 
tures. Only when the forward scattering character of the 
impurity potential is further enhanced, the vertex correc- 
tions begin to outweigh the growing self-energy and the 
conductivity rises above the values obtained for isotropic 
scattering. For these more extended scattering potentials 
the second resonance in the self-energy moves to lower 
frequencies until it merges with first resonance. For the 
larger microwave frequency of 13.4 GHz, see right panel 
of Fig. the anisotropy of the impurity potential has 
less effect. This implies that for slightly extended im- 
purity potentials as considered for YBCO in Fig. [3] the 
frequency dependence becomes much weaker and there- 
fore a smaller concentration than in the case of pointlikc 
impurities is necessary to reproduce the experimentally 
observed frequency dependence. We emphasize that this 
analysis cannot rule out other explanations^ 2 ^ of the 
quasilinear in T behavior observed at low frequencies. 



IV. COMPARISON WITH EXPERIMENTAL 
SPECTRA OF BSCCO 

In this section we want to explore what can be learned 
about the type of disorder contained in the BSCCO com- 
pounds by analyzing its microwave conductivity which 
has measured by S.-F. Lee et alH. The temperature and 
frequency dependence of the microwave conductivity in 
BSCCO (see symbols in Fig. is quite different than in 
YBCO. The absolute value of the microwave conductivity 
is smaller by almost a factor of 10 indicating that BSCCO 
is a dirtier compound than YBCO. This agrees well with 
the observation that the conductivity of BSCCO changes 
noticeably only in the THz-regimeA^, i.e., at much higher 
frequencies than in YBCO. The characteristic peak in the 
microwave conductivity, however, occurs at lower tem- 
peratures than in the cleaner system YBCO contrary to 
predictions for strong pointlike scatterers 2 . Furthermore 
this peak is much less pronounced and resembles more 
a plateau, suggestive of the weak scattering limi1j£. Fi- 
nally, the conductivity does not apparently approach the 
universal value ao for the lowest temperatures and fre- 
quencies measured. This might indicate the presence of 
extended scatterers, which enhance the zero-temperature 
and zero-frequency limit of the conductivity as has been 
shown by Durst and Lee22. 

In order to check the applicability of these scenar- 
ios for BSCCO we compare in Fig. [7| respective fits 
to the microwave conductivity using (i) only pointlike 
weak scatterers (Fig. Ufa)) and (ii) only extended weak 
scatterers (Fig. Efb)). For the inelastic scattering rate 
T inci(T) we assume the same temperature-dependence as 
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FIG. 7: Fit to the experimental microwave conductivity of 
BSCCO (reproduced from Ref. ll2fl using weak scatterers: (a) 
pointlike scatterers with V/t = 1 and n = .049, (b) weak 
extended scatterers with V\/t=2, V2/t=.8, Vz/t=A and n — 
.145. (The inelastic scattering rate has been increased by a 
factor of 2.6 in (a) and 3 in (b) with respect to Ref . I29T) . 



for YBCOS but we allow for a different prefactor in order 
to account for discrepancies between YBCO and BSCCO 
(the exact prefactor is stated in the figure captions). Ob- 
viously both disorder models (i) and (ii) result in very 
good fits of the experimental microwave conductivity of 
BSCCO. The main difference is that the conductivity for 
isotropic scatterers approaches the universal value do for 
T — > whereas the conductivity for extended scatter- 
ers (Fig. EJb)) with the potential parameters V\/t = 2, 
Vz/t = 0.8 and V^/t — 0.4 approaches an enhanced value. 
Unfortunately, it is not possible to distinguish between 
these two different scenarios from the microwave conduc- 
tivity alone, because there is no experimental data avail- 
able below T = 5K. 

Further insight could be gained by comparing the fre- 
quency dependence of these two models. Impurity scat- 
tering alone would predict a different frequency depen- 
dence for isotropic and extended scatterers, as illustrated 
in Fig.|Hl Whereas the magnitude of the conductivity re- 
mains rather large at low temperatures in the case of 
extended scatterers even for high frequencies, it almost 
vanishes in the case of pointlike impurities. Due to the 
large impurity concentration, however, this frequency de- 
pendence is most pronounced in the THz-regime, as ob- 
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FIG. 8: Higher frequency conductivity for the disorder models 
of Fig. [T] (a) pointlike scatterers with V/t = 1 and n = .049, 
(b) weak extended scatterers with Vi/t=2, l^/t=.8, Vz/t=A 
and n = .145. 



served experimentally^, where the contribution of inelas- 
tic scattering cannot be neglected. Thus, the THz-data 
present not only a probe of elastic impurity scattering 
but also of inelastic scattering processes and are there- 
fore not directly suitable to distinguish between pointlikc 
and extended impurities. 

The only way we can proceed now is try to exclude one 
of the two models indirectly via analyzing an additional 
observable. Thus we will argue in the following that dis- 
order model (i) containing 4.9% weak isotropic scatterers 
with a scattering strength of V=lt would yield an unre- 
alistically large normal state scattering rate. Assuming 
a normal state density of states of 1 /St yields a normal 
state scattering rate of t _1 ps 0.6T c for t=120meV. Ac- 
cording to Abrikosov Gorkov's scaling law this would 
reduce T c by 25%, which we consider as an unreason- 
ably large suppression because T c w 93K in the samples 
used for microwave conductivity in Ref. ll2T which is close 
to the highest values of T c measured for the BSCCO- 
compounds. Extended impurities, on the other hand, 
act mainly as small angle scatterers and affect T c much 
less than isotopic scatterers^*^. This allows us to assess 
model (i), which consists only of weak pointlike scatter- 
ers, as very unlikely and to conclude that at least a large 
fraction of the disorder in the BSCCO compounds should 
be attributed to extended scatterers. This could be con- 
firmed by experiments on crystals at lower T. 

So far we have focused on the effect of weak scatterers 
which we attribute to the out-of plane disorder intro- 
duced by doping. Defects whiting the Cu02-planes like 
Zn-substitution for Cu, on the other hand, are generally 
considered to act as strong pointlike scatterers. These 
strong scattering defects have been observed in STM- 
experiments^ on BSCCO-compounds and are often as- 
sumed to be the main source of disorder in the YBCO- 
compounds although their concentration is very low. We 
therefore now address the question of whether our model 
for the microwave conductivity in BSCCO is compatible 
with an additional small concentration of strong point- 
like impurities, which are most likely also present in the 
compound used for measurements of the microwave con- 
ductivity of BSCCO in Ref.[T3 

We incorporate this realistic disorder model, which 
consists of weak extended and strong pointlike scatter- 
ers, by calculating the diagrams depicted in Fig. |5J The 
weak extended scatterers are treated in Born approxi- 
mation and the strong pointlike impurities in T-matrix 
approximation, where as before we consider only the So- 
component of the self-energy, i.e., only diagrams of the 
T-matrix with an even number of impurity lines. Because 
vertex corrections vanish for pointlike scatterers, only the 
weak extended scatterers contribute and the vertex cor- 
rections can be calculated in Born approximation. Thus, 
the microwave conductivity cr(f2) is still given by Eq. 1L'2I) 
with the bubble J{lo, £1) as in Eq. (|18|) an d the vertex 
function j BoTn given in Eq. I|23|) . Only the self-energy 
Eq has to be calculated self-consistently as the sum of 
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FIG. 9: Diagrams for the self-energy E and microwave con- 
ductivity a considered in the realistic disorder model for 
BSCCO. Circles denote the pointlike strong impurities, which 
are treated in T-matrix approximation. Crosses denote the 
weak extended scatterers, which are treated in Born approx- 
imation. Only the extended, i.e., only the weak scatterers 
contribute to the vertex corrections. 
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Eq. i|12|) . Eq. (|14|) and the inelastic scattering rate 

Eo(w < T) = -Ii^- i-wW? ) (25) 

-n w {V? + 2V 2 2 + V 3 2 )I g (lj) - z(2r ine i(T))- 1 . 

Here, n s is the concentration of strong pointlike impu- 
rities with a scattering potential V s , n w is the concen- 
tration of weak extended scatterers characterized by the 
potential parameters V\, V%, V3 and Ig(u) is the momen- 
tum integrated single particle Green's function given by 
Eq. ©. 

The effect of adding a small concentration of strong 
pointlike scatterers to the extended weak scatterers used 
in Fig. EJb) is illustrated in Fig. ^| Additional strong 
pointlike scatterers mainly reduce the conductivity as can 
be seen Fig. HUf a) . In order to raise the conductivity to 
its previous values the forward scattering character of 
the weak impurities must therefore be enhanced. On the 
other hand, this increases the difference between the zero- 
temperature value and the maximum value of the conduc- 
tivity (see dashed line in Fig. HOf b)) and necessitates a 
larger concentration of weak extended scatterers. Finally 



FIG. 11: Fit to the experimental micowave conductivity of 
BSCCO (reproduced from Ref. 112) using n w = .1 weak ex- 
tended scatterers with Vi/t—3, V2/t=.9, Vs/t=.5 and n s = 
.0005 strong pointlike scatterers with V s /t = 100. (The in- 
elastic scattering rate has been increased by a factor of 3.4 
with respect to Ref. I29I) 



the lineshape of the conductivity (solid line in Fig. llOf b)) 
looks similar to the one without strong pointlikc impu- 
rities but it is flatter at low temperatures than before. 
This poses an upper limit for the concentration of strong 
scatterers compatible with the experimental microwave 
conductivity of BSCCO. 

A fit to the experimental microwave conductivity of 
BSCCO containing .05% pointlike strong scatterers and 
10% weak extended scatterers is shown in Fig.^] This is 
about the largest concentration of pointlike strong scat- 
terers still compatible with the microwave conductivity. 
This concentration is smaller than the .2% observed in 
STM-experiments 23 on BSCCO but it is very plausible 
that the number of in-plane defects varies for different 
compounds and maybe even between bulk and surface. 
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FIG. 10: Effect of adding strong pointlike scatterers for the 
microwave conductivity at 14.4 GHz. (a) Weak extended 
scatterers with parameters of Fig. Efb), i.e., n w — 0.145 and 
Vi/t = 2, Vilt = 0.8, Vzjt = 0.4, and additional strong point- 
like scatterers with V/t — 100 and different concentrations 
n s . (b) Fixed concentration n s = 0.001 of strong pointlike 
impurities for varying concentration and forward scattering 
potential of the weak scatterers. 



V. CONCLUSIONS 

In summary, we have investigated the influence of 
extended scatterers on the microwave conductivity of 
d-wave superconductors by extending the approach of 
Durst and Lee^S, which is based on a nodal approxi- 
mation for the quasiparticle spectrum and the impurity 
potential, to finite temperatures and frequencies. We 
have slightly modified the approach of Durst and Lee by 
considering only vertex functions which contain an even 
number of impuritity lines. This modification is neces- 
sary to ensure analyticity at finite temperatures and fre- 
quencies when only the normal self-energy component So 
is considered. 

The effect of extended scatterers on the temperature 
and frequency dependence of the microwave conductiv- 
ity can be summarized as follows: for a small concentra- 
tion of slightly extended strong scatterers a second reso- 
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nance forms in the self-energy at intermediate frequencies 
similar to treatments which consider the suppression of 
the order parameter surrounding a strong pointlike impu- 
rity^. This results in a more linear temperature depen- 
dence of the conductivity at low temperatures and there- 
fore improves the agreement with experimental spectra of 
YBCO at low temperatures. For more extended scatter- 
ers the vertex corrections begin to dominate over the self- 
energy and the magnitude of the conductivity increases. 

The microwave conductivity of BSCCO is very differ- 
ent compared to YBCO and cannot be understood in 
terms of only strong scattering pointlikc impurites. We 
find that a large concentration of weak extended scat- 
terers is necessary to explain the observed temperature 
and frequency dependence of the microwave conductiv- 
ity in BSCCO, where: (i) the impurity concentration 
has to be large to explain the small magnitude of the 
conductivity and the negligible frequency dependence in 
the microwave range, (ii) the scattering strength has to 
be small to account for the plateau-like lineshape of the 
conductvity at small temperatures and (iii) the range of 



the scattering potential has to be spatially extended in 
order to keep the T c -suppression reasonably small^i. Fi- 
nally, we have shown that adding a small concentration of 
pointlike strong scatterers, which have been observed in 
STM-experiments^, to the weak extended scatterering 
potential, which we attribute to the out-of-plane disor- 
der introduced by doping, is still compatible with the 
microwave conductivity of BSCCO. Although it would 
be necessary to refine our treatment of inelastic scatter- 
ing by accounting for its frequency dependence and its 
contribution to vertex corrections in order to address the 
conductivity in the THz-rangei^, it is interesting to note 
that elastic scattering alone would predict a very dif- 
ferent frequency dependence for pointlike and extended 
scatterers. 
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